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Abstract. We survey the theory of increasing and decreasing subsequences of permuta¬ 
tions. Enumeration problems in this area are closely related to the RSK algorithm. The 
asymptotic behavior of the expected value of the length is{w) of the longest increasing 
subsequence of a permutation w of 1, 2,... ,n was obtained by Vershik-Kerov and (al¬ 
most) by Logan-Shepp. The entire limiting distribution of is(u;) was then determined by 
Baik, Deift, and Johansson. These techniques can be applied to other classes of permuta¬ 
tions, such as involutions, and are related to the distribution of eigenvalues of elements of 
the classical groups. A number of generalizations and variations of increasing/decreasing 
subsequences are discussed, including the theory of pattern avoidance, unimodal and 
alternating subsequences, and crossings and nestings of matchings and set partitions. 


1. Introduction 

Let &n denote the symmetric group of all permutations of [n] := {1, 2, ..., n}. We write 
permutations w £ &„ as words, i.e., w = 0102 • • • dn, where ui(i) = at. An increasing 
subsequence of w is a subsequence ■ ■ ■ Oi^, satisfying < • • • < ai ^,, and similarly for 
decreasing subsequence. For instance, if w = 5642713, then 567 is an increasing subse¬ 
quence and 543 is a decreasing subsequence. Let is(u;) (respectively, ds(u;)) denote the 
length (number of terms) of the longest increasing (respectively, decreasing) subsequence 
of w. If rc = 5642713 as above, then is(w) = 3 (corresponding to 567) and ds(rc) = 4 (cor¬ 
responding to 5421 or 6421). A nice interpretation of increasing subsequences in terms of 
the one-person card game patience sorting is given by Aldous and Diaconis [3|. Further 
work on patience sorting was undertaken by Burstein and Lankham EH 113 Connec¬ 
tions between patience sorting and airplane boarding times were found by Bachmat et 
al. I11II12I and between patience sorting and disk scheduling by Bachmat 1101 . 

The subject of increasing and decreasing subsequences began in 1935, and there has 
been much recent activity. There have been major breakthroughs in understanding the 
distribution of is(u;), ds(r(;), and related statistics on permutations, and many unexpected 
and deep connections have been obtained with such areas as representation theory and 
random matrix theory. A number of excellent survey papers have been written on various 
aspects of these developments, e.g., El EH] EH Gnu mil; the present paper will emphasize 
the connections with combinatorics. 

In Section 1^ we give some basic enumerative results related to increasing/decreasing 
subsequences and show their connection with the RSK algorithm from algebraic com¬ 
binatorics. The next two sections are devoted to the distribution of is(u;) for w £ &n, 
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a problem first raised by Ulam. In Section we deal with the expectation of is(w) for 
ui € 6n, culminating in the asymptotic formula of Logan-Shepp and Vershik-Kerov. We 
turn to the entire limiting distribution of is(w) in Section 0] The main result is the de¬ 
termination of this limiting distribution by Baik, Deift, and Johansson to be a (suitably 
scaled) Tracy-Widom distribution. The Tracy-Widom distribution originally arose in the 
theory of random matrices, so the result of Baik et al. opens up unexpected connections 
between increasing/decreasing subsequences and random matrices. 

Much of the theory of increasing/decreasing subsequences of permutations in 
carries over to permutations in certain subsets of &n, such as the set of involutions. 
This topic is discussed in Section|^ In particular, analogues of the Baik-Deift-Johansson 
theorem were given by Baik and Rains. In Section|^we explain how the previous results 
are related to the distribution of eigenvalues in matrices belonging to the classical groups. 

The remaining three sections are concerned with analogues and extensions of increas¬ 
ing/decreasing subsequences of permutations. Section |3 deals with pattern avoidance, 
where increasing/decreasing subsequence are replaced with other “patterns.” In Section|^ 
we consider unimodal and alternating subsequences of permutations, and in Section|^we 
replace permutations with (complete) matchings. For matchings the role of increasing 
and decreasing subsequences is played by crossings and nestings. 

Acknowledgment. I am grateful to Percy Deift, Persi Diaconis, Craig Tracy and 
Herb Wilf for providing some pertinent references. 


2. Enumeration and the RSK algorithm 

The hrst result on increasing and decreasing subsequences is a famous theorem of Erdos 
and Szekeres mi 

Theorem 1. Let p, g > 1. If w £ ©p^+i, then either is(w) > p or ds(w) > g. 

This result arose in the context of the problem of determining the least integer f(n) 
so that any f(n) points in general position in the plane contain an n-element subset S 
in convex position (i.e., every element of S' is a vertex of the convex hull of S). A recent 
survey of this problem was given by Morris and Soltan m Seidenberg gave an 
exceptionally elegant proof of Theorem 0 based on the pigeonhole principle which has 
been reproduced many times, e.g., Gardner Ml Ch. 11, §7]. 

Theorem^is best possible in that there exists w £ &pq with is{w) — p and ds(u)) = g. 
Schensted IHHI found a quantitative strengthening of this result based on his rediscovery 
of an algorithm of Robinson m which has subsequently become a central algorithm in 
algebraic combinatorics. To describe Schensted’s result, let A = (Ai, A 2 ,...) be a partition 
of n > 0, denoted A h n or |A| = n. Thus Ai > A 2 > • • • > 0 and Xi = n. The (Young) 
diagram of a partition A is a left-justified array of squares with Ai squares in the ith row. 
For instance, the Young diagram of (3, 2, 2) is given by 
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A standard Young tableau (SYT) of shape A h n is obtained by placing the integers 
1,2, ...,n (each appearing once) into the squares of the diagram of A (with one integer 
in each square) such that every row and column is increasing. For example, an SYT of 
shape (3, 2, 2) is given by 

136 

25 . 

47 

Let denote the number of SYT of shape A. The quantity has a number of 
additional combinatorial interpretations Ell Prop. 7.10.3]. It also has a fundamental al¬ 
gebraic interpretation which suggests (via equation below) a close connection between 
representation theory and increasing/decreasing subsequences. Namely, the (complex) ir¬ 
reducible representations of 6n are indexed in a natural way by partitions A h n, and 
then 

f^ = dimF\ (1) 

In particular (by elementary representation theory), 

E (h)’ = »!. (2) 

Ahn 

See Section]^ for more on the connections between increasing/decreasing subsequences 
and representation theory. MacMahon El P- 175] was the first to give a formula for 
(in the guise of counting “lattice permutations” rather than SYT), and later Frame, 
Robinson, and Thrall 1431 simplified MacMahon’s formula, as follows. Let u be a square 
of the diagram of A, denoted u £ A. The hook length h(u) of (or at) u is the number of 
squares directly to the right or directly below u, counting u itself once. For instance, if 
A = (3, 2, 2) then the hook lengths are given by 

541 

32 . 

21 

The hook-length formula of Frame, Robinson, and Thrall asserts that if A h n, then 


For instance. 


j(3,2,2) 


7! 

5-4-1-3-2-2-1 


= 21 . 


The RSK algorithm gives a bijection between ©„ and pairs (P, Q) of SYT of the 
same shape A h n. This algorithm is named after Gilbert de Beauregard Robinson, who 
described it in a rather vague form El §5] (subsequently analyzed by van Leeuwen El 
§7]), Craige Schensted 1361 . and Donald Knuth El For further historical information 
see El Notes to Ch. 7]. The basic operation of the RSK algorithm is row insertion, i.e., 
inserting an integer i into a tableau T with distinct entries and with increasing rows and 
columns. (Thus T satisfies the conditions of an SYT except that its entries can be any 
distinct integers, not just 1,2, ...,n.) The process of row inserting i into T produces 
another tableau, denoted T ^ i, with increasing rows and columns. If S is the set of 
entries of T, then S U {i} is the set of entries oi T i. We define T *— i recursively as 
follows. 




4 


• If the first row of T is empty or the largest entry of the first row of T is less than 
i, then insert i at the end of the first row. 

• Otherwise, i replaces (or bumps) the smallest element j in the first row satisfying 
j > i. We then insert j into the the second row of T by the same procedure. 

For further details concerning the definition and basic properties of T <— f see e.g. EH 
Ch. SJEl §7.11]. 

Let w = aia 2 • • • a„ G 6„, and let 0 denote the empty tableau. Define 
Pi = Pi{w) = (■ • • ((0 ^ ai) ^ 02 ) ^ ai. 

That is, we start with the empty tableau and successively row insert 01 , 02 ,... ,ai. Set 
P — P{w) — Pn{w). Define Qo = 0, and once Qi-i is defined let Qi = Qi{'w) be obtained 
from Qi-i by inserting i (without changing the position of any of the entries of Qi-i) 
so that Qi and Pi have the same shape. Set Q = Q{w) = Qniw), and finally define the 
output of the RSK algorithm applied to w to be the pair {P,Q), denoted w ^ {P,Q). 
For instance, if w = 31542 G 65 , then we obtain 

Pi (in) = 1, P 2 {w) = J , P 3 (w) = 3 ® 

14 12 
Pa{w) ^ 35 . P ^ Ps{w) = 3 4 . 

5 

It follows that 

13 

Q= 24 . 

5 

Note. By a theorem of Schutzenberger EHEH §7.13] we have 

Q{w) = P(w“^), (4) 

so we could have in fact taken this formula as the definition of Q{w). 

If w ^ (P, Q) and P, Q have shape A, then we also call A the shape of w, denoted 
A = sh(«;). The conjugate partition A' = (Ai, A 2 ,...) of A is the partition whose diagram 
is the transpose of the diagram of A. Equivalently, j occurs exactly Xj — Aj+i times 
as a part of A'. The length i{X) is the number of (nonzero) parts of A, so ^(A) = A). 
The fundamental result of Schensted EH connecting RSK with increasing and decreasing 
subsequences is the following. 

Theorem 2. Let w G &n, and suppose that sh(w) = A. Then 

is(u)) = Ai (5) 

ds(ui) = A 1 . ( 6 ) 

Equation 0 is easy to prove by induction since we need only analyze the effect of the 
RSK algorithm on the first row of the Pi’s. On the other hand, equation 0 is based on 
the following symmetry property of RSK proved by Schensted. If w = 0102 • • • a„ then 
let = On - ■ ■ 0201 , the reverse of w. We then have 

{P,Q) ^ w’'^ (P‘,evac((3)‘), 


(7) 
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where * denotes transpose and evac{Q) is a certain tableau called the evacuation of Q 
(first defined by Schfitzenberger ED) which we will not define here. Equation 0 shows 
that if sh(r(;) = A, then sh(w’^) = A'. Since clearly is{w^) = ds(M;), equation 0 follows 
from 0 . 

Theorem |21 has several immediate consequences. The first corollary is the Erdos- 
Szekeres theorem (Theorem 0, for if sh(w) = A, is(ui) < p, and ds(u;) < q, then Ai < p 
and Ai < q. Thus the diagram of A is contained in a, q x p rectangle, so |A| < pq. By the 
same token we get a quantitative statement that Theorem 0 is best possible. 

Corollary 3. The number of permutations w G Gpq, where say p < q, satisfying is(u;) = 
p and ds(u;) — q is given by 

/r(p «)'|2 ^ ( _(m)]_ 

\ 1^22 • • . pP(p + l)p . • • qP(q + l)p-i ■ • • (p + g — 1)1 

where (p®) denotes the partition with q parts equal to p. 

Proof. Let A = sh(M;). If is(ui) = p and ds(M;) = g, then Ai = p and A'l = g. Since 
A h pq, we must have A = (p'*). The number of u € 6 n with a fixed shape p is just 
the number of pairs (P, Q) of SYT of shape p. Hence the left-hand side of equation 0 
follows. The right-hand side is then a consequence of the hook-length formula 0 . □ 

An interesting result concerning the extremal permutations in the case p = g in 
Corollary 0 was obtained by Romik IS2I Thm. 5]. It can be stated informally as follows. 
Pick a random permutation w £ ©^2 satisfying is(M;) = ds(M;) = p. Let Pw be the p^ x p^ 
permutation matrix corresponding to w, drawn in the plane so that its corners occupy 
the points (±1,±1). Then almost surely as p —> oo the limiting curve enclosing most of 
the I’s in P™ is given by 

{(i, p) G : (x^ - y^f + 2{x^ -I- p^) = 3}. 

See Figure 0 In particular, this curve encloses a fraction a = 0.94545962 • • • of the 
entire square with vertices (±1, ±1). The number a can be expressed in terms of elliptic 
integrals as 



a = 2 


/■ 


v/(l-t2)(i-(f/3)2) 


zdt — 


-T 

2 Jo 


1 - {t/^y 


dt. 


Compare with the case of any w G &n, when clearly the limiting curve encloses the entire 
square with vertices (±1,±1). For further information related to permutations w G 6^2 
satisfying is(w) = ds(w) = p, see the paper EH of Pittel and Romik. 

Clearly Corollary0can be extended to give a formula Cor. 7.23.18] for the number 


Qpqin) = #{w G 6 „ : is(u;) = p, ds(ui) = g}. 


namely, 

9p<iin) = ■ ( 9 ) 

Ahn 

>■ 1 =?, >>'1=9 

The usefulness of this formula may not be readily apparent, but Theorem 0 below is an 
example of its utility. 
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3. Expectation of is{w) 

A further application of Theorem concerns the distribution of the function is(w) as 
w ranges over The problem of obtaining information on this distribution was first 
raised by Ulam [THSl §11.3] in the context of Monte Carlo calculations. In particular, one 
can ask for information on the expectation E{n) of is(w) for w £ i.e., 

Ulam mentions the computations of E. Neighbor suggesting that E{n) is about 1.7^/n. 
Numerical experiments by Baer and Brock m suggested that E(n) ~ 2yTi might be 
closer to the truth. The Erdos-Szekeres theorem (Theorem^ implies immediately that 
E{n) > y/n, since 

i(is(ui) +is(ui’')) > ^/is{w)\s{nF) = \/is(u)) ds(w) > V^- 

Hammersley [EH] was the first person to seriously consider the question of estimating 
E{n). He showed that if 

hm 

n—*-oo 

then c exists and satisfies 

TT 

2 <c<e. 

He also gave a heuristic argument that c = 2, in agreement with the experiments of Baer 
and Brock. 
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The next progress on Ulam’s problem was based on Schensted’s theorem (Theorem|^. 
It follows from this result that 


= (/')'• ( 10 ) 

Ahn 

Now the RSK algorithm itself shows that 

= ( 11 ) 

Ahn 

in agreement with Since the number of terms in the sum on the right-hand side of 
itTTl is very small compared to n!, the maximum value of for A h n is close to -hnl. 
Let A" be a value of A h n for which is maximized. Then by TO we see that a close 
approximation to E{n) is given by 

«(A'*)i. 


This heuristic argument shows the importance of determining the partition A" maximizing 
the value of for Ahn. 

We are only really interested in the behavior of A" as n —> oo, so let us normalize 
the Young diagram of any partition A to have area one. Thus each square of the diagram 
has length Xjypn. Let the upper boundary of (the diagram of) A be the y-axis directed 
to the right, and the left boundary be the x-axis directed down. As n —> oo it is not 
unreasonable to expect that the boundary of the partition A"' will approach some limiting 
curve y = Vl'(a;). If this curve intersects the z-axis at a: = 6 , then it is immediate that 


c := 



> b. 


We cannot be sure that 6 = c since conceivably the first few parts of A'* are much larger 
than the other parts, so these parts would “stretch out” the curve y = along the 

a;-axis. 

It was shown independently by Vershik-Kerov nna and Logan-Shepp 1681 that y = 
^'(a:) indeed does exist and is given parametrically by 


X = y + 2 cos 9 
2 

y = — (sin0 — 9 cos 6), 

TV 

ioT 0 < 9 < TV. See Figure]^ where we have placed the coordinate axes in their customary 
locations. 

Logan-Shepp and Vershik-Kerov obtain the curve (a;) as a solution to a variational 
problem. If {x,y) is a point in the region A enclosed by the curve and the coordinate 
axes, then the (normalized) hook-length at {x,y) is f{x) — y + f~^{y) — x. By equation 
0 we maximize by minimizing Hug a ^(w). Hence in the limit we wish to minimize 
I (/) = //^ log(/(a;) — y + f~^{y) — x)dx dy, subject to the normalization fdx dy — 1. 
It is shown in m and Enzi that / = is the unique function minimizing /(/) (and 
moreover /(/) = — 1 / 2 ). 
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1.5 

1 

0.5 

0 



Figure 2. The curve y = (a:) 


Note. If we extend the curve y — T(a:) to include the a;-axis for x > 2 and the t/-axis 
for y > 2, and then rotate it 45° counterclockwise, then we obtain the curve 


Q{x) = 


|■( 2 ; arcsin(a;/2) + Vd — |a:| < 2, 

l»l, l®l > 2. 


This form of the limiting curve is more convenient for some purposes EUEOIEI]. such as 
surprising connections with the separation of zeros of orthogonal polynomials. 

We see immediately from the equations for T( 2 :) that it intersects the i-axis at a: = 2, 
so c > 2. By a simple but clever use of the RSK algorithm Vershik and Kerov show in 
their paper that c < 2, so we conclude that 


E{n) ~ 2\/n. 


( 12 ) 


Different proofs that E{n) ~ 2-^/n were later given by Aldous and Diaconis 0, Groene- 
boom m, Johansson and Sepalainen IM]. The proof of Aldous and Diaconis is 
based on an interacting particle process for which the number of particles remaining af¬ 
ter n steps has the same distribution as iSn. Their proof is known in the language of 
statistical physics as a hydrodynamic limit argument. See |11 §3] for a brief survey. 

We should remark that the curve y = 4'(a:) is not merely the limiting curve for the 
partition maximizing /^; it is also the limiting curve for the typical shape of a permuta¬ 
tion w € &n- A remarkable refinement of this fact is due to Kerov I53|[5H][SI1 §0.3.4], who 
shows that the deviation of a Young diagram from the expected limit converges in prob¬ 
ability to a certain Gaussian process. A different kind of refinement is due to Borodin, 
Okounkov, and Olshanski EHI Thm. 1], who obtain more detailed local information about 
a typical shape A than is given by T(a;). 
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4. Distribution of is{w) 

A major breakthrough in understanding the behavior of is{w) was achieved in 1999 by 
Baik, Deift, and Johansson El They determined the entire limiting distribution of is(u;) 
as n —> 00 . It turns out to be given by the (suitably scaled) Tracy-Widom distribution, 
which had first appeared in connection with the distribution of the largest eigenvalue of 
a random hermitian matrix. 

To describe these results, write is„ for the function is : ©„ —> Z. Let u{x) denote the 
unique solution to the nonlinear second order differential equation 

u" {x) = 2u{x)^ + xu{x), (13) 


subject to the condition 


u(x) 






00 . 


Equation GSJ is known as the Painleve II equation, after Paul Painleve (1863-1933). 
Painleve completely classified differential equations (from a certain class of second or¬ 
der equations) whose “bad” singularities (branch points and essential singularities) were 
independent of the initial conditions. Most of the equations in this class were already 
known, but a few were new, including equation GSJ. 

Now define the Tracy-Widom distribution to be the probability distribution on R 
given by 


H 


F{t) = exp - / (a; — t)u(x) dx 


(14) 


It is easily seen that F{t) is indeed a probability distribution, i.e., F{t) > 0 and F'{t)dt 
1. We can now state the remarkable results of Baik, Deift, and Johansson. 


Theorem 4. We have for random (uniform) w £ &n and all t £ R that 


lim Prob 

n —>00 


/ iSn(rc) — 2y/n 



= m- 


The above theorem is a vast refinement of the Vershik-Kerov and Logan-Shepp results 
concerning E{n), the expectation of is(M;). It gives the entire limiting distribution (as 
n —> oo) of isTi(w). Baik, Deift, and Johansson also determine all the limiting moments 
of is„('u;). In particular, we have the following formula for the variance Var(is„) of is„ as 
n —» CX3. 


Corollary 5. We have 


and 


= 0.8131947928 • • • , 

lim - f tdF{t) (15) 

n—*co 71^'^ J 

= -1.7710868074-•• . 
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Note that equation may be rewritten 

E{n) = 2^/n + an}^^ + o(n^'^®), 

where a = / tdF{t), thereby giving the second term in the asymptotic behavior of E{n). 

We will say only a brief word on the proof of Theorem|l| explaining how combinatorics 
enters into the picture. Some kind of analytic expression is needed for the distribution 
of iSn{w). Such an expression is provided by the following result of Ira Gessel 1451 . later 
proved in other ways by various persons; see ITU §1] for references. Define 


Uk{n) = #{w e 6 „ : is„(w) < k} 

'^k{n)—pr, k>l 
n\^ 

n>0 


h{2x) = ^ 

n>0 


n\ {n + i)\’ 


i e Z. 


(16) 


The function E is the hyperbolic Bessel function of the first kind of order i. Note that 
Ii{2x) = I-i{2x). 

Theorem 6. We have 

Uk{x) = det {Ii-j{2x))^^ .^^. 

Example 1. We have (using /i = I-i) 


U2{x) = 


Io{2.x) 
h{2x) 

= h{2xf- 


h{2x) 

Io{2x) 

Ii{2xf. 


From this expression it is easy to deduce that 


M 2 (n) 


1 

n + 1 



(17) 


a Catalan number. This formula for U 2 {n) was first stated by Hammersley 1581 in 1972, 
with the first published proofs by Knuth 1031 §5.1.4] and Rotem m There is a more 
complicated expression for U 3 (n) due to Gessel EH §7)1001 Exer. 7.16(e)], namely. 


U3{n) 


1 ^ /2 A /n + 1 Wn + 2\ 

(n + l)2(n + 2) jyj + ljyj + 2j' 


(18) 


while no “nice” formula for Uk{n) is known for fixed fc > 3. It is known, however, that 
Uk{n) is a P-recursive function of n, i.e., satisfies a linear recurrence with polynomial 
coefficients EOl §7]. For instance, 

(n + 4)(n + 3 )^U 4 (n) = (20n® + 62n^ + 22n — 24)u4(n — 1) 

— 64n(n — l)^U 4 (n — 2) 


(n + 6)^(n + Afu^in) = (375 - 400n - 843n^ - 322n^ - 35n'‘)u5(n - 1) 
+ (259n^ + 622n + 45)(n - ifusin - 2) 

- 225(n - l)^(n - 2fu5{n - 3). 
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A number of conjectures about the form of the recurrence satisfied by Uk{n) were made 
by Bergeron, Favreau, and Krob m, reformulated in with some progress toward a 
proof. 

Gessel’s theorem (Theorem reduces the theorem of Baik, Deift, and Johansson 
to “just” analysis, viz., the Riemann-Hilbert problem in the theory of integrable sys¬ 
tems, followed by the method of steepest descent to analyze the asymptotic behavior of 
integrable systems. For further information see the survey m of Deift mentioned above. 

The asymptotic behavior of is„ (w) (suitably scaled) turned out to be identical to the 
Tracy-Widom distribution F{t) of equation 11411 . Originally the Tracy-Widom distribu¬ 
tion arose in connection with the Gaussian Unitary Ensemble (GUE). GUE is a certain 
natural probability density on the space of all n x n hermitian matrices M = (Mij), 
namely, 

where is a normalization constant and 

dM = n dMii ■ d(Re Mij)d{lm Mij). 

i i<j 

Let the eigenvalues of M be ai > a 2 > The following result marked the 

eponymous appearance UMI of the Tracy-Widom distribution: 

lim Prob — \/^') \/2n^^® < t) — F{t). (19) 

Thus as n —> oo, iSn{w) and a\ have the same distribution (after scaling). 

It is natural to ask, firstly, whether there is a result analogous to equation II 9> for 
the other eigenvalues ak of the GUE matrix M, and, secondly, whether there is some 
connection between such a result and the behavior of increasing subsequences of random 
permutations. A generalization of cu to all Qfc was given by Tracy and Widom umi 
(expressed in terms of the Painleve II function u(x)). The connection with increasing 
subsequences was conjectured in m and proved independently by Borodin-Okounkov- 
Olshanski isni, Johannson m, and Okounkov ES], after first being proved for the second 
largest eigenvalue by Baik, Deift, and Johansson ca Given w G &n, define integers 
Ai, A 2 ,... by letting Ai -I- ■ • • -I- A^ be the largest number of elements in the union of k 
increasing subsequences of w. For instance, let w = 247951368. The longest increasing 
subsequence is 24568, so Ai = 5. The largest union of two increasing subsequences is 
24791368 (the union of 2479 and 1368), so A 1 + A 2 = 8 . (Note that it is impossible to find a 
union of length 8 of two increasing subsequences that contains an increasing subsequence 
of length Ai = 5.) Finally w itself is the union of the three increasing subsequences 
2479, 1368, and 5, so Ai -I- A 2 -I- A 3 = 9. Hence (Ai, A 2 ,A 3 ) = (5,3,1) (and Ai = 0 for 
i > 3). Readers familiar with the theory of the RSK algorithm will recognize the sequence 
(Ai, A 2 ,...) as the shape sh(ui) as defined preceding Theorem El a well-known result of 
Gurtis Greene 14811961 Thm. Al.1.1]. (In particular, Ai > A 2 > • • • , a fact which is by 
no means obvious.) The result of 1 : lull 'll iim asserts that as as n ^ cxa, Afc and Ok are 
equidistributed, up to scaling. In particular, the paper of Okounkov provides a direct 
connection, via the topology of random surfaces, between the two seemingly unrelated 
appearances of the Tracy-Widom distribution in the theories of random matrices and 
increasing subsequences. A very brief explanation of this connection is the following: 
a surface can be described either by gluing together polygons along their edges or by 
a ramified covering of a sphere. The former description is related to random matrices 
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via the theory of quantum gravity, while the latter can be formulated in terms of the 
combinatorics of permutations. 

We have discussed how Gessel’s generating function Uk{x) for Uk{n) is needed to find 
the limiting distribution of is„. We can also ask about the behavior of Uk(n) itself for 
fixed k. The main result here is due to Regev m 


Theorem 7. For fixed k and for n —> oo we have the asymptotic formula 


Uk{n) ~ 1! 2! • • • (fc — 1)! 




,y/2 


k^ 


Idea of proof. From the RSK algorithm we have 


,(fc2-l)/2- 


Uk 


(")= E (/“)^ 


( 20 ) 


Ahn 


Write in terms of the hook-length formula, factor out the dominant term from the 
sum (which can be determined via Stirling’s formula), and interpret what remains in the 
limit n —» 00 as a fc-dimensional integral. This integral turns out to be a special case of 
Selberg’s integral (e.g., Ch. 8]), which can be explicitly evaluated. □ 

An immediate corollary of Theorem [7| (which can also be easily proved directly using 
RSK) is the formula 

lim Mfc(n)^^'‘ = k^. (21) 


5. Symmetry. 

Previous sections dealt with properties of general permutations in Much of the 

theory carries over for certain classes of permutations. There is a natural action of the 
dihedral group Di of order 8 on ©n, best understood by considering the permutation 
matrix corresponding to w G &n- Since is a square matrix, acts on P™ as the 
usual symmetry group of the square. In particular, reflecting through the main diagonal 
transforms P-w to its transpose P,^ = Pw-^- Reflecting about a horizontal line produces 
Pw’’' 5 v^here ru is the reverse of w as used in equation 0. These two reflections generate 
the entire group D/^. 

Let G be a subgroup of Di, and let 

6 ® = {w € &n : cr • ui = w for all cr € G}. 

Most of the results of the preceding sections can be carried over from 6„ to ©)^. The 
general theory is due to Baik and Rains [DiniiiH] Moreover, for certain G we can add 
the condition that no entry of P^, equal to 1 can be fixed by G, or more strongly we can 
specify the number of I’s in P^, fixed by G. For instance, if G is the group of order 2 
generated by reflection through the main diagonal, then we are specifying the number of 
fixed points of w. For convenience we will consider here only two special cases, viz., (a) 
G is the group of order 2 generated by reflection through the main diagonal. In this case 
= {w € ©n : w'^ = 1}, the set of involutions in ©„, which we also denote as 3n. (b) 
The modification of (a) where we consider fixed-point free involutions only. Write 3^ for 
this set, so 3)) = 0 when n is odd. 
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The RSK algorithm is well-behaved with respect to inversion, viz., it follows from 
equation 0 that if ui ^ {P,Q) then w~^ ^ {Q,P). Hence 

ui^ = 1 if and only if P = Q. (22) 


Let 


yk{n) = #{w € : is„(w) < fc}. 

By Schensted’s theorem (Theorem |2J we conclude 


yk{n) = 

\\-n 

Xi<k 


the “involution analogue” of JUJ. From this formula or by other means one can ob¬ 
tain formulas for yk{n) for small k analogous to itTTt and itTSl . In particular (see 
Exer. 7.16(b)] for references), 

2/2 (n) 

2/3 (n) 

2 / 4 (^^) 

2/5(41) 

where as usual Ci is a Catalan number. 

The RSK algorithm is also well-behaved with respect to fixed points of involutions. 
It was first shown by Schiitzenberger 1871 p. 1271 1961 Exer. 7.28(a)] that if = 1 and 
w (P, P), then the number of fixed points of w is equal to the number of columns of 
P of odd length. Let 



« 2 fc(n) = #{w £ Dn : ds(w) < 2k} 

Zk{n) = e : is(«;) < k}. 

(It is easy to see directly that if w £ then ds(M;) is even, so there is no need to deal 
with V 2 k+iin).) It follows that 


V2k{n) = 

Ahn 

Ai<fc 

Zk{n) = Y ’ 

Ahn 
A^ <fe 

where 2A' = (2Ai, 2 A 2 ,...), the general partition with no columns of odd length. Note 
that for fixed-point free involutions w £ we no longer have a symmetry between is(rc) 
and ds(M;), as we do for arbitrary permutations or arbitrary involutions. 
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There are also “involution analogues” of Gessel’s determinant (Theorem]^. Equations 
and {HI below were first obtained by Gessel SSI §6], equation ll^ by Goulden m, 
and equations {HI and m by Baik and Rains |1 fil Gor. 5.5]. Let 

n 

Yk{x) = — 

^^ n! 

n>0 

V2kin)^ 

ni 

n>0 

Zhin) — . 
n\ 

n>0 


Write li = Ii{2x). 

Theorem 8. We have 


y^ 2 k{^) — det {li—j + — 

Y 2 k+i{x) = det {li-j - 

V 2 k{x) = det {li—j — — l 

Z 2 k{x) = i det {h-j + + i det {h-j - 

Z 2 k+i{x) = ie^ det {h-j - det {h-j + Ii+j-i)\ 




(23) 

(24) 

(25) 

(26) 
(27) 


Once we have the formulas of Theorem]^ we can use the techniques of Baik, Deift, 
and Johansson to obtain the limiting behavior of ds(u)) for ui G and w € These 
results were first obtained by Baik and Rains cum 

Theorem 9. (a) We have for random (uniform) w & 'Jn and alH G R that 


lim Prob 

n—>oo 


is„(w) — 2y/n 

ni /6 


< t ) = F(t)^/^exp 


1 r 

2 A 


u{s)ds , 


where F{t) denotes the Tracy-Widom distribution and u{s) the Painleve II function. (By 
{HI we can replace isn(w) with dsn(w).) 

(b) We have for random (uniform) w £ ^271 and all t G R that 


lim Prob (< t) = exp 

n^oo (2n)i/® - ' ^ ' 


1 

2Jt 


u(s)ds . 


(c) We have for random (uniform) w G J 2 n and all t G R that 


lim Prob 

n —>00 


/ iS2n(w) — 2V^ 
(2n)i/6 


< t ) = E(t)^/^cosh 


1 r 
2Jt 


u{s)ds j . 


There are orthogonal and symplectic analogues of the GUE model of random hermi- 
tian matrices, known as the GOE and GSE models. The GOE model replaces hermitian 
matrices with real symmetric matrices, while the GSE model concerns hermitian self-dual 
matrices. (A 2n x 2n complex matrix is hermitian self-dual if is composed of 2 x 2 blocks of 
a bi c-\- di 


the form 


—c -i- di a — bi 


which we identify with the quaternion a-\-bi-\-cj -\-dk, such 
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that if we regard the matrix as an n x n matrix M of quaternions, then Mji = Mij where 
the bar is quaternion conjugation.) The limiting distribution of ds 2 n{w) for w G 
coincides (after scaling) with the distribution of the largest eigenvalue of a random real 
symmetric matrix (under the GOE model), while the limiting distribution of is 2 n(rc) 
for w G ^271 coincides (after scaling) with the distribution of the largest eigenvalue of a 
random hermitian self-dual matrix (under the GSE model) 11041 . 


6. Connections with the classical groups. 


In equation 0 we expressed gpq(n), the number of permutations w G ©n satisfying 
is(rc) = p and ds(u;) = q, in terms of the degrees of irreducible representations of 6„. 
This result can be restated via Schur-Weyl duality as a statement about the distribution 
of eigenvalues of matrices in the unitary group U{n). The results of Section|^can be used 
to extend this statement to other classical groups. 

Let U{k) denote the group of k x k complex unitary matrices. For a function / : 
U{k) —> C, let E{f) denote expectation with respect to Haar measure, i.e., 

E{f) = f f{M)dM, 

J MeU(k) 

where f is the Haar integral. The following result was proved by Diaconis and Shahsha- 
hani m for n > k and by Rains m for general k. Note that if M has eigenvalues 
01 ,...,6k then 

|tr(M)"|= = (6)i + ... + 6k)’^{0q + ... + 4)T 

Theorem 10. We have i5(|tr(M)"|^) = Uk{n), where Uk{n) is defined in equation 11611 . 

Proof. The proof is based on the theory of symmetric functions, as developed e.g. in m 
or EHl Gh. 7]. If fixi ,..., Xk) is a symmetric function, then write /(M) for f{9i ,..., 6k), 
where 9i,... ,9k are the eigenvalues of M G U{k). The Schur functions sx for ^(A) < k 
are the irreducible characters of U{k), so by the orthogonality of characters we have for 
partitions A, p of length at most k that 

/ sx{M)7;jJf)dM = Sxp. (28) 

J MeU(k) 

Now tr(M)" = pi(M)", where pi(a:i,..., Xk) = xi + ■ ■ ■ + Xk. The symmetric function 
Pi{xi,... ,Xk)" has the expansion 1691 Exam. 1.5.21 1??^ Cor. 7.12.5] 


Piixi 


,Xk 


)" = X] f^sxixi,...,Xk), 


Ahn 

i{X)<k 


where is the number of SYT of shape A as in Section |21 (This formula is best under¬ 
stood algebraically as a consequence of the Schur-Weyl duality between ©„ and U{k) 

Ch. 7, Appendix 2], although it can be proved without any recourse to representation 
theory.) Hence from equation 1281 we obtain 

E(|tr(M)"|^) = [ pi(M)>i(M)"dM 
J MGU{k) 

- E {f'T- 

Ahn 

^(A)<fc 
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Comparing with equation completes the proof. □ 

Many variations of Theorem llOlhave been investigated. For instance, we can replace 
tr(M)" by more general symmetric functions of the eigenvalues, such as tr(M"*)", or 
we can replace U{k) with other classical groups, i.e., 0{k) and Sp(2fc). For further 
information, see Rains ca 


7. Pattern avoidance. 


In this and the following two sections we consider some generalizations of increasing/de¬ 
creasing subsequences of permutations. In this section and the next we look at other 
kinds of subsequences of permutations, while in Section|^we generalize the permutations 
themselves. 

We have defined Uk{n) to be the number of permutations in &n with no increasing 
subsequence of length k + 1. We can instead prohibit other types of subsequences of a 
fixed length, leading to the currently very active area of pattern avoidance. 

Given v = bi ■ ■ ■ bk G &k, we say that a permutation w = ai ■ ■ ■ On € &n avoids v if 
it contains no subsequence • ai^. in the same relative order as u, i.e., no subsequence 
ciii • ■ • Oif, satisfies the condition: 

V 1 < r < s < fc, br < bs- 


Thus a permutation w satisfies is(w) < fc if and only if it is 12 ■ • • fc-avoiding, and similarly 
satisfies ds(ui) < fc if and only if it is fc(fc — 1) • • ■ 1-avoiding. What can be said about the 
set &n{v) of permutations w € &n that are u-avoiding? In particular, when are there 
formulas and recurrences for Sn{v) := #&n{v) similar to those to those of Theorem]^ 
and Example 0? 

The vast subject of pattern avoidance, as a generalization of avoiding long increasing 
and decreasing subsequences, began in 1968 with Knuth 1681 Exer. 2.2.1.5]. He showed 
in connection with a problem on stack sorting that Sn(312) is the Catalan number Cn- 
(See also Exer. 6.19(ff)].) By obvious symmetries this result, together with equation 
o. shows that Sn{v) = Cn. for all v € 63. A fundamental paper directly connecting 
321-avoiding and 231-avoiding permutations was written by Simion and Schmidt m- 

Wilf first raised the question of investigating Sn(u) for v & &k when fc > 4. Here is 
a brief summary of some highlights in this burgeoning area. Eor further information, see 
e.g. Hi Chs. 4,5] and the special issue [H|. Call two permutations u,v £ &k equivalent, 
denoted u ~ w, if Sn{u) = Sn{v) for all n. Then there are exactly three equivalence 
classes of permutations u £ 64. One class contains 1234, 1243, and 2143 (and their 
trivial symmetries), the second contains 3142 and 1342, and the third 1324. The values 
of Sn(1234) are given by equation 11811 and of Sn(1342) are given by the generating function 


^Sn(1342)a;" 

n>0 


_ 32 ®_ 

1 + 20 ® — 8 ®^ — (1 — 8 ®)®/^ ’ 


a result of Bona m .The enumeration of 1324-avoiding permutations in ©„ remains open. 

Let us mention one useful technique for showing the equivalence of permutations in 
6 fc, the method of generating trees introduced by Chung, Graham, Hoggatt and Kleiman 
m and further developed by West irn^irnnuTTiTi and others. Given u £ Gk, the gen¬ 
erating tree Tu is the tree with vertex set IJn>i and with y a descendent of w if 
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4312 3412 3142 4132 1432 4321 3421 

4312 3412 3124 4123 1234 4321 3421 


3241 3214 4231 2431 4213 2413 2143 

3241 3214 2413 2134 4231 2341 2314 


Figure 3. The generating tree for ©„(123) and 6„(132) 


u) is a subsequence of y (an actual subsequence, not the pattern of a subsequence). For 
many pairs w,n G 6fc we have Tu =%, showing in particular that Sn{u) = Sn(v) for all 
n, i.e., u ~ V. In many cases in fact the two trees will have no automorphisms, so the 
isomorphism —> 77 is unique, yielding a canonical bijection 6n{u) —> &n{v). A unique 

isomorphism holds for instance when u = 123 and v = 132. Figure|^shows the first four 
levels of the trees 7i23 = 7i32, labelled by elements of both 7i23 and 7i32 (boldface). This 
tree can also be defined recursively by the condition that the root has two children, and 
if vertex x has k children, then the children of x have 2, 3,..., fc + 1 children. For further 
information about trees defined in a similar recursive manner, see Banderier et al. m 
Given v G &k, let 

Fv{x) = ^ s„(v)x'‘. 

n>0 

It is not known whether F7(a:) is always algebraic or even the weaker condition of being 
D-finite, which is equivalent to Sn{v) being P-recursive 19311961 §6.4]. A long-standing 
conjecture, known as the Stanley-Wilf conjecture, stated that for all v G &k there is a 
c > 1 such that Sn{v) < c". In this case lim„^oo Sn(v)^^" exists and satisfies 

1 < < 00 0. For instance, equation (HU cLSSBftS tllclt Zjl2- -k = (fc — 1)^. The Stanley- 

Wilf conjecture was proved by Marcus and Tardos by a surprisingly simple argument. 
It is known that for fc > 4 the permutation 12 - ■ ■ k neither maximizes nor minimizes 
for i> € ©fc [2, but it is not known which permutations do achieve the maximum or 
minimum. 

An interesting aspect of pattern avoidance was considered by Albert Let X he a, 
finite subset of ©2 U ©3 U ■ • ■ , and let Sn{X) denote the number of permutations w G &n 
avoiding all permutations v G X. We say that X is proper if it doesn’t contain both the 
identity permutation 12 • ■ • j for some j and the “reverse identity” k - ■ - 21 for some k. It 
is easy to see (using Theorem 0 that X is proper if and only if Sn{X) > 0 for all n > 1. 
For w G ©n let is^(u;) be the length of the longest subsequence of w avoiding all v G X, 
and let Ex{ri) denote the expectation of is^(w) for uniform w G ©n, so 

Ex{n) = is^(«’)- 

w£&n 

Albert then makes the following intriguing conjecture. 
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Conjecture 11. If X is proper then 


lim i 

n—*oo 4 



Ex{n) \ 


2 


(29) 


(The limit on the left-hand side exists by m and a simple generalization of [7|, while 
the limit on the right-hand side is only conjectured to exist.) 

Consider for instance the case X = {21}. Then Sti(-Y) = 1 for all n > 1, so the 
left-hand side of is 1. On the other hand, is^(w) = is(w), so the right-hand side is 
also 1 by equation Gl. More generally, Albert proves Conjecture mi for X = {12---k} 
and hence (by symmetry) for X = {k - ■ ■ 21} ^ Prop. 4]. 

Let us mention that permutations avoiding a certain pattern u or a finite set of 
patterns have arisen naturally in a variety of contexts. For instance, an elementary 
result of Tenner ]m asserts that the interval [0, ui] in the Bruhat order of &n is a 
boolean algebra if and only if w is 321 and 3412-avoiding, and that the number of such 
permutations w is the Fibonacci number Fbn-i. The Schubert polynomial is a single 
monomial if and only if w is 132-avoiding COl p. 46]. All reduced decompositions of 
w G &n are connected by the Coxeter relations SiSj = SjSi (i.e., SiSi+iSi does not appear 
as a factor in any reduced decomposition of w) if and only if w is 321-avoiding 1241 
Thm. 2.1]. Vexillary permutations may be defined as those permutations w such that the 
stable Schubert polynomial FL is a single Schur function or equivalently, whose Schubert 
polynomial is a flag Schur function (or multi-Schur function). They turn out to be the 
same as 2143-avoiding permutations ITni (1.27)(iii), (7.24)(iii)], first enumerated by West 
ESHl Cor. 3. 1711117111 Cor. 3.11]. Similarly, permutations w G Sn for which the Schubert 
variety in the complete flag variety GL(n, C)/B is smooth are those permutations that 
are 4231 and 3412-avoiding (implicit in Ryan 1841 . based on earlier work of Lakshmibai, 
Seshadri, and Deodhar, and explicit in Lakshmibai and Sandhya EB]). The enumeration 
of such “smooth permutations” in is due to Haiman 12811501117(71 Exer. 6.47], viz.. 


^ 6n(4231,3412)a;”' 

n>0 


1 — X 


1 


( _2x_ 

1 — x y 1+a; —(1 —a:)C(a:) 



where C{x) = Yln>o = (1 — \/l — 4:x)/2x. See Billey and Lakshmibai 1251 for 

further information. As a final more complicated example, Billey and Warrington m 
show that a permutation w € &n has a number of nice properties related to the Kazhdan- 
Lusztig polynomials Px,ni if and only if w avoids 321, 46718235, 46781235, 56718234, and 
56781234. These permutations were later enumerated by Stankova and West 1921 . A 
database of “natural occurrences” of pattern avoidance can be found at a website m\ 
maintained by B. Tenner. 

The subjects of pattern avoidance and increasing/decreasing subsequences can be 
considered together, by asking for the distribution of is(w) or ds(ui) where w ranges over 
a pattern-avoiding class Gn(v)- (For that matter, one can look at the distribution of is(w) 
or ds(rc) where w ranges over any “interesting” subset of 0^.) For instance, Reifegerste 
uni Cor. 4.3] shows that for k > 3, 


#{m; e 6 to( 231) : is(M;) < fc} 




(30) 


a sum of Narayana numbers 1961 Exer. 6.36]. Note that the left-hand side of 131)1 can also 
be written as #{w £ ©n(231,12 • ■ • fc)}, the number of permutations in ©n avoiding both 
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231 and 12 - ■ ■ k. Asymptotic results were obtained by Deutsch, Hildebrand, and Will m 
for the distribution of is(ui) when v = 231,132, and 321. Their result for v = 132 is the 
following. 

Theorem 12. For w G G„(132) the random variable is(u;) has mean y/nn + 
and standard deviation — l)\/u + Moreover, for any t > —y/n we have 

lim Prob ~ ^ < t] =y(l-2/(f + 0F)2)e-(‘+'^)'^'. 

n-i-oo Y -Vn 7 

\ y / HczT 


The proof of Theorem ll!2l is considerably easier than its counterpart for w G &n (Theo¬ 
rem because there is a relatively simple formula for the number f{n, k) of permutations 
w G 6n(132) satisfying is{w) < k, viz., 


f{n,k) = 2 


L(n+l)/(fc+l)J 

E 

I=r-n/(fe+l)l 


2 n 

n + i(k + 1) 


1 / 2n + 2 \ 

4 (n + 1 + i(fc + l)y 


A number of variations and generalizations of pattern-avoiding permutations have 
been investigated. In particular, we can look at patterns where some of the terms must 
appear consecutively. This concept was introduced by Babson and Steingrimsson and 
further investigated by Claesson and others. For instance, the generalized pattern 
1-32 indicates a subsequence aiajaj+i of a permutation w = aia 2 ■ ■ ■ an such that ai < 
aj+i < aj. The hyphen in the notation 1-32 means that the first two terms of the 
subsequence need not be consecutive. The permutations in G 4 avoiding 1-32 are all 
C4 — 14 permutations avoiding 132 (in the previous sense, so avoiding 1-3-2 in the 
present context) together with 2413. A typical result, due to Claesson 1371 Props. 2 and 
5], asserts that 

#e„(l-23) = #e„(l-32) = B{n), 
the number of partitions of the set [n] (a Bell number ESI §l-4])- 


8. Unimodal and alternating subsequences. 

We briefly discuss two variations of increasing/decreasing subsequences of a different 
flavor from those considered above. There is considerable room for further work in this 
area. 

Early work of Chung ES] and Steele m deals with fc-unimodal subsequences. A 
sequence is k-unimodal if it is a concatenation of (at most) k + 1 monotone sequences. (In 
more traditional terminology, a sequence is fc-unimodal if it has at most fc -|- 1 alternating 
runs EHl §i- 21 .) Thus a 0-unimodal sequence is just an increasing or decreasing sequence, 
and every such sequence is fc-unimodal for all k. The sequences 41235 and 24531 are 1- 
unimodal. Chung showed that every w £ &n has a 1-unimodal subsequence of length 
-k — ij , and that this result is best possible. She conjectured that if Ek{n) is 
the expected length of the longest fc-unimodal subsequence of a random permutation 
w £ Gn, then Ek{n)/y/n approaches a limit Ck as n ^ 00 . Steele proved this conjecture 
and showed that Ck = 2y/k + 1 by deducing it from the monotone {k = 0 ) case. 

A sequence 6162 • • • 6 fc of integers is alternating if 


bi > b2 < bs > bi < ■ ■ ■ bk- 
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For instance, there are five alternating permutations in © 4 , viz., 2143, 3142, 4132, 3241, 
4231. If En denotes the number of alternating permutations in &n, then a famous result 
of Andre ElESl §3.16] states that 


n>0 


X 

n\ 


secx + tana:. 


(31) 


The numbers En were first considered by Euler (using m as their definition) and are 
known as Euler numbers. Sometimes E 2 n is called a secant number and E 2 n-i a tangent 
number. 

We can try to extend the main results on increasing/decreasing subsequences to alter¬ 
nating subsequences. In particular, given w £ ©^ let as{w) = asn(rc) denote the length 
of the longest alternating subsequence of w, and define 


bk{n) = #{w £ ©„ : as(w) < k}. 

Thus bi (n) = 1 (corresponding to the permutation 12 • • • n), bk (n) = n! if fc > n, and 
bn{n) — bn-i(n) = En. Note that we can also define bk{n) in terms of pattern avoidance, 
viz., bk{n) is the number of re £ ©„ avoiding all Ek+i alternating permutations in ©fc+i. 

Unlike the situation for Uk{n) (defined by llfit i. there are “nice” explicit generating 
functions and formulas for bk{n). The basic reason for the existence of such explicit 
results is the following (easily proved) key lemma. 


Lemma 13. For any w € Gn, there exists an alternating subsequence of w of maximum 
length that contains n. 

From Lemma ESI it is straightforward to derive a recurrence satisfied by ak{n) := 
bk{n) - bk-i{n), viz., 

= E ( 1 ) E ia2rij-1)+a2r+iij-l))asin-j). (32) 

j = l / 2r + s = fc-l 


Now define 

n 

B{x,t)= bk{n)t'" — . 

nl 

k,n>0 

It follows from the recurrence 13211 (after some work IhHI I that 


B{x, t) 


1 + p + 2te'’=^ + {1 - p)e^'’^ 

1 + p — t^ + {l — p — t^)e^P^ ’ 


(33) 


where p = y/1 — B. Alternatively (as pointed out by M. Bona), let G{n,k) denote the 
number of w £ &n with k alternating runs as defined at the beginning of this section. 
Then equation 1331 is a consequence of the relation ak(n) = ^{G{n, fc — 1) -|- G(n, k)) and 
known facts about G{n, k) summarized in El §l- 2 ]. 

It can be deduced from equation 13311 (shown with assistance from I. Gessel) that 


bk{n) 


2k-l 


E (-2^ 


i+2j<k 
= k (mod 2) 


k- j 
{k + x )/2 


For instance, 

&i(n) = l, & 2 (n)= 2 "-\ fe3(n) = i(3"- 


2n-|-3), fo 4 (n) = ■i(4"'— (2n — 4)2"). 

O 
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From equation (El it is also easy to compute the moments 

Mfc(n) = -7 as„(«;)'". 

r) I • ^ 


For instance, 


^ Mi(r 


dB{x, t) 


dt 


6x — 3x^ + 

6(1 — x)'^ 


from which we obtain 


Mi{n) ^ I 4 n + l 


6 


n = 1 


n > 1 . 


Similarly the variance of asTi('U’) is given by 
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var(as„) = —n-n > 4. 

^ ^ 45 180 - 


(34) 


(35) 


It is surprising that there are such simple explicit formulas, in contrast to the situation 
for is(w) (equation lllOll l. 

It is natural to ask for the limiting distribution of asn, analogous to TheoremEJfor isn- 
The following result was shown independently by R. Pemantle IZl and H. Widom ma¬ 
lt can also be obtained by showing that the polynomials ak(n)t^ have (interlacing) 
real zeros, a consequence of the connection between ak(n) and G{n, k) mentioned above 
and a result of Wilf El Thm. 1.41], 

Theorem 14. We have for random (uniform) w € &n and all t G R that 

lim p^^^'as^(w)-2n/3 
n —>00 


< t = G(t), 


where G{t) is Gaussian with variance 8/45: 




tV45/4 


ds. 


For further information on longest alternating subsequences, see the paper 


9. Matchings. 

The subject of pattern containment and avoidance discussed in Section [3 provides one 
means to extend the concept of increasing/decreasing subsequences of permutations. In 
this section we will consider a different approach, in which permutations are replaced with 
other combinatorial objects. We will be concerned mainly with (complete) matchings on 
[2n], which may be defined as partitions M = {Bi,..., B^} of [2n] into n two-element 
blocks Bi. Thus Bi U 52 U ■ • • U Bn = [2n], Bi n Bj = 0 if i / j, and #Bi = 2. 
(These conditions are not all independent.) Alternatively, we can regard a matching M 
as a fixed-point free involution wm of [2n], viz., if Bi = {a,b} then WM{a) = h. We 
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already considered increasing and decreasing subsequences of fixed-point free involutions 
in Section 13 In that situation, however, there is no symmetry interchanging increasing 
subsequences with decreasing subsequences. Here we consider two alternative statistics 
on matchings (one of which is equivalent to decreasing subsequences) which have the 
desired symmetry. 

Write Tin for the set of matchings on [2n]. We represent a matching M € Tin by a 
diagram of 2n vertices 1, 2 ,..., 2n on a horizontal line in the plane, with an arc between 
vertices i and j and lying above the vertices if {i,j} is a block of M. Figure^Jshows the 
diagram corresponding to the matching 

M = {{1, 5}, {2, 9}, {3,10}, {4, 8}, {6, 7}}. 

Let M e Tin- A crossing of M consists of two arcs {i, j} and {k, 1} with i < k < j < 1. 
Similarly a nesting of M consists of two arcs {i,j} and {k,l} with i < k < I < j. The 
maximum number of mutually crossing arcs of M is called the crossing number of M, 
denoted cr(M). Similarly the nesting number ne(M) is the maximum number of mutually 
nesting arcs. For the matching M of Figure 2] we have cr(M) = 3 (corresponding to the 
arcs {1,5}, {2,9}), and {3, 10}), while also ne(M) = 3 (corresponding to {2,9}, {4,8}, 
and {6, 7}). 

It is easy to see that ds(uiM) = 2 • ne(M), where wm is the fixed-point free involution 
corresponding to M as defined above. However, it is not so clear whether cr(M) is 
connected with increasing/decreasing subsequences. To this end, define 

Mi,j) = #{M € Tin ■■ cr(M) = i, ne(M) = j}. 

It is well-known that 

= (36) 

j » 

In other words, the number of matchings M £ Tin with no crossings (or with no nest¬ 
ings) is the Catalan number Cn- For crossings this result goes back to Errera 14211961 
Exer. 6.19(n,o)]; for nestings see Ea Equation 11361 was given the following generalization 
by Chen et al. 1341 . 

Theorem 15. For all i,j,n we have fn{i,j) = fn{j,i)- 

Theorem IT3 is proved by using a version of RSK first defined by the author (un¬ 
published) and then extended by Sundaram 111)61 . Define an oscillating tableau of shape 
A h n and length fc to be a sequence 

0 = A°,A\...,A'' =A 
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Figure 5. 

A correspondence between 

matchings and oscillating tableaux 


of partitions A* such that (the diagram of) A'^^ is obtained from A® by either adding or 
removing a square. (Note that if we add a square each time, so k = n, then we obtain an 
SYT of shape A.) Oscillating tableaux were first defined (though not with that name) by 
Berele QH in connection with the representation theory of the symplectic group. Given 
a matching M G DTln, define an oscillating tableau ${M) = (A°, A^,..., A^") of length 2n 
and shape 0 as follows. Label the right-hand endpoints of the arcs of M by 1, 2,..., n 
from right-to-left. Label each left-hand endpoint with the same label as the right-hand 
endpoint. Begin with the empty tableau To = 0. Let ai,..., a 2 n be the sequence of labels, 
from left-to-right. Once Ti_i has been obtained, define Ti to be the tableau obtained by 
row-inserting ai into Ti-\ (as defined in Section|2l if is the label of a left-hand endpoint 
of an arc; otherwise T is the tableau obtained by deleting ai from T-i. Let A® be the 
shape of T, and set 

$(M) = (0 = A°,A\...,A^" =0). 

See Figure for an example. It is easy to see that ^{M) is an oscillating tableau of 
length 2n and shape 0. With a little more work it can be shown that in fact the map 
M <1>(M) is a bijection from Tin to the set On of all oscillating tableaux of length 2n 
and shape 0. As a consequence we have the enumerative formula 

#On = (2n - 1)!! ;= 1 ■ 3 • 5 • • • (2n - 1), (37) 

the number of matchings on [2n]. The key fact about the correspondence $ for proving 
Theorem 1 151 is the following “oscillating analogue” of Schensted’s theorem (Theorem I^J. 

Theorem 16. Let M G Tin, and suppose that <I>(M) = (A°, A^,..., A^"). Then ne(M) is 
equal to the most number of columns of any A\ while cr{M) is equal to the most number 
of rows of any A*. 

Theorem 1151 is an easy consequence of Theorem 1161 For let ^{M)' be the oscillating 
tableau obtained by conjugating all the partitions in $(M), and let M' be the matching 
satisfying ${M') = <1?(M)'. By Theorem 1161 we have cr(M) = ne(M') and ne(M) = 
cr(M'). Since the map M i—> <1?(M) is a bijection we have that the operation M i—> M' is 
a bijection from Tin to itself that interchanges cr and ne, and the proof follows. 

The above argument and equation Q show that the operation M i—> M' on matchings 
is a natural analogue of the operation of reversal on permutations. Unlike the case of 
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permutations, we don’t know a simple “direct” operation on matchings that interchanges 
cr with ne. 

Let fn denote the number of oscillating tableaux of shape A and length n, so fn = 0 
unless n = |A| (mod2). A generalization for of the hook-length formula (equation l|^) 
is due to Sundaram [MD Lemma 2.2], viz.. 


P = 

J n 





where A h fc and where of course is evaluated by the usual hook-length formula 
(Set (—1)!! = 1 when n = k.) Note that an oscillating tableau (A°, A^,..., A^") of shape 
0 may be regarded as a pair (P, Q) of oscillating tableaux of the same shape A = A" and 
length n, viz., 

P=(A°,A\...,A") 

Q = (A2",a2'*-\...,A"). 

Hence we obtain the following restatement of equation P7|l: 

^(/„^)' = (2n-l)!!, (38) 


where A ranges over all partitions. (The partitions A indexing a nonzero summand are 
those satisfying A h fc < n and k = n (mod 2).) 

Equation I38II suggests, in analogy to equations Q and Q, a connection between 
and representation theory. Indeed, there is a C-algebra 23^(a;), where a; is a real 
parameter, which is semisimple for all but finitely many x (and such that these exceptional 
X are all integers) and which has a basis that is indexed in a natural way by matchings 
M e In particular, dim®TO(a;) = (2n — 1)!!. This algebra was first defined by Brauer 
m and shown to be the centralizer algebra of the action of the orthogonal group 0(H) 
on H®" (the nth tensor power of V), where dim V — k and x = k. It is also the centralizer 
algebra of the action of the symplectic group Sp(2fc) on H®", where now dimH = 2k 
and X = —2k. When 23^(a:) is semisimple, its irreducible representations have dimension 
fn, so we obtain a representation-theoretic explanation of equation iT^ . For further 
information see e.g. Barcelo and Ram 1201 App. B6] 

Because ds(u;M) = 2 • ne(M) and because cr„ and nen have the same distribution by 
Theorem uni the asymptotic distribution of cr„ and nCn on reduces to that of ds 2 n 
on 32n, which is given by Theorem |^b). We therefore obtain the following result. 

Theorem 17. We have for random (uniform) M G Tin and all t £ R that 


lim Prob 

n —>-00 


/nen(M) — ^ 

(2n)i/6 - 2 ) 



The same result holds with ne„ replaced with cr„. 


We can also consider the effect of bounding both cr(M) and ne(M) as n —» cx3. The 
analogous problem for is(w) and ds(ui) is not interesting, since by Theorem Q there are 
no permutations w £ &n satisfying is(ui) < p and ds(u;) < q as soon as n > pq. Let 


hp,q{n) = #{M £ Mn ■ cr(M) < p, ne{M) < q} 
Hp,q(x) = ^/ip,9(n)a:". 

n>0 
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It follows from the bijection $ : DJtn —> On and a simple application of the transfer-matrix 
method §4.7] that Hp^q{x) is a rational function of x §5]. For instance, 




l-2x 
1 — 3a: -I- 


„ , . 1 — 5z -I- 2x^ rr / \ 1 “ 11® -I- 30a;^ — 23a:® -f 4x'* 

il-a:){l-5xy (l-:r)(l-3a:)(l-8a: + 4a:2) 

1 - 243: -f 186x^ - 567a:® -f 690a:'^ - 285x® -f ISa:*^ 

(1 — a:)(l — 19a: -f 83x^ — x^){l — 5x -I- 6x^ — a:^)^ 

Christian Krattenthaler pointed out that hp^qin) can be interpreted as counting certain 
walks in an alcove of the affine Weyl group Cn- It then follows from a result of Grabiner 
m (23)1 that all reciprocal zeros of the denominator of Hp^q{x) are of the form 


2(cos(7rri/m) -I- • • • -I- cos{'Krj/m)), 


where each G Z and m = p + q + 1. All these numbers for fixed m belong to an extension 
of Q of degree <^(2m)l2, where 0 is the Euler phi-function. As a consequence, every 
irreducible factor (over Q) of the denominator of Hp^q{x) has degree dividing (f>{2m)/2. 

Theorem llhi can be extended to objects other than matchings, in particular, arbitrary 
set partitions. (Recall that we have dehned a matching to be a partition of [2n] into n 2- 
element blocks.) In this situation oscillating tableaux are replaced by certain sequences of 
partitions known as vacillating tableaux. See 1341 for further details. Vacillating tableaux 
were introduced implicitly (e.g., [521 (2.23)]) in connection with the representation theory 
of the partition algebra a semisimple algebra whose dimension is the Bell number 
B{n). See Halverson and Ram 1521 for a survey of the partition algebra. Vacillating 
tableaux and their combinatorial properties were made more explicit by Chen, et al. m 
and by Halverson and Lewandowski [51]. An alternative approach based on “growth 
diagrams” to vacillating tableaux and their nesting and matchings was given by Krat¬ 
tenthaler [55]. It remains open to End an analogue of Theorem 1 171 for the distribution of 
cr(7r) or ne(7r) (as defined in [M]) for arbitrary set partitions tt. 
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